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We investigate the short-time dynamics of a delta-function potential barrier on an initially confined 
wave-packet. There are mainly two conclusions: A) At short times the probability density of the 
first particles that passed through the barrier is unaffected by it. B) When the barrier is absorptive 
(i.e., its potential is imaginary) it affects the transmitted wave function at shorter times than a real 
potential barrier. Therefore, it is possible to distinguish between an imaginary and a real potential 
barrier by measuring its effect at short times only on the transmitting wavefunction. 

PACS numbers: 03.65.-w, 03.65.Nk, 03.65.Xp. 
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I. INTRODUCTION 



One of the methods in imaging a turbid or a diffusive 
medium with optical radiation is the time-gating tech- 
nique PISH SOI- I n this technique a temporally narrow 
pulse is injected into the medium. Owing to the diffu- 
sivity of the medium, when the pulse exits the medium 
it becomes considerably wider. However, if the first ar- 
riving photons are separated from the rest of the pulse 
then it is possible to use these, so called, ballistic (or 
quasi ballistic) photons to reconstruct the ballistic image 
of the medium. 

Therefore, by employing a short time-gating tech- 
nique the multi-scattering effect can be eliminated. In- 
deed, such methods were employed in recognizing hid- 
den objects and informative signals in diffusive media 
Naively, one might expect that this technique can 
be implemented for electron imaging to "see" absorp- 
tive objects inside scattering medium. That is, one can 
send a short pulse of electrons to one end of the medium, 
while at the other end only the first arriving electrons 
will be measured. By doing so all the noise caused by 
the multi-scattering should be eliminated. On the other 
hand, the presence of absorptive regions (like imaginary 
potentials 0) will be felt in the amount of the early 
arriving electrons. However, electrons are governed by 
the Schrodinger equation, and unlike the Maxwell wave- 
equation, has a parabolic dispersion relation. As a con- 
sequence, any localized wave-packet suffers from strong 
dispersion, since each spectral component propagates at 
different velocity. The fastest particles are the most en- 
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ergetic ones, which pass through the medium unaffected, 
since the barrier's potential energy is negligible compared 
to their kinetic energy. In particular, when the medium 
is a certain barrier (or well), we conclude that when only 
the first-arriving particles are measured there should be 
no trace of the barrier's presence. It does not matter 
what shape or height the barrier has - the first parti- 
cles that pass though the barrier should be indifferent to 
it. Thus, one may argue that the time-gating technique 
cannot be implemented to electrons imaging, at least not 
in its naive form. However, we show that the short-time 
measurement reveals information about the nature of the 
barrier- whether it is imaginary or real. 

There is a peculiar distinction between an absorb- 
ing medium (e.g., an imaginary potential) and a non- 
absorbing one (e.g., a real potential). While they both 
have no effect on the wavepacket (both transmitted and 
reflected) at t — > 0, the imaginary barrier influences the 
wavepacket sooner. In other words, in the temporal Tay- 
lor expansion of the probability density the imaginary 
potential appears at smaller order than a real potential. 
It is then clear that we can classify the barrier as an 
absorptive one simply by measuring the wave-packet at 
short times. Note that in general it is required to measure 
both reflection and transmission coefficients to figure out 
if the barrier is absorptive or not. 

Recently |8| , it has been demonstrated even experimen- 
tally that it is feasible to investigate the ID scattering of 
a Bose-Einstein condensate by a narrow defect. There- 
fore, it seems that there is a good chance of witnessing 
these effects in the laboratory in the near future. 

In this paper we demonstrate this effect rigorously 
(both analytically and by a numerical simulation) for the 
delta function potential. That is, we show that it is pos- 
sible to identify an absorptive potential by measuring the 
short time dynamics of only the transmitted wavefunc- 
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FIG. 1: System schematic: 
a delta function potential. 



tion. 



a semi-infinite wavepacket hitting 



The initial state we consider is a wave-packet, which 
is confined to one side of the barrier. It is then demon- 
strated that the wavefunction at the other side is inde- 
pendent of the barrier for short times, while the temporal 
dependence depends on the exact nature of the barrier 
(absorptive or not). 



II. SYSTEM DESCRIPTION AND DYNAMICS 

Evidently, in order to confine the initial wavepacket to 
one side of the barrier, there has to be a certain singular- 
ity in the wavepacket. In this paper we focus on a step 
function to simplify the problem; however, it has been 
demonstrated elsewhere that most of the conclusions are 
valid even in the continuous case provided the measure- 
ments are taken at specific ranges (see ref. 9j). Is is also 



demonstrated at the end of the paper that the main con- 
clusions are valid even when the initial wavepacket is a 
Gaussian. For simplicity we take a ID delta function as 
the potential barrier. 

The system illustration is depicted in fig.l. Initially, 
the wave packet has the form 0, 0, [H| 

ip(x, t = 0) = 9 (— x) exp (ikox) (1) 

and a distance L from its front we place a delta- 
function barrier V(x) — XS(x — L) (see fig.l). That is, 
the Schrodingcr equation reads 
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therefore for t > 
ip(x,t) 
where 

<p(k) 



dkcp(k)x(k 1 x)exp(-ik 2 t) (3) 
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k — ko + i0 
is the Fourier transform of ip(x, t = 0) and 



X(k,x) = exp(ikx) + - 



iX/2 



k-iX/2 

The general solution is (for x > 0) 



exp (ik\x — L\ + ikL) (5) 
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where w(z) = exp(— z 2 )erfc(— iz) [TlT ] and y = L + \x — L\. In short times one can expand this expression in powers 
of t. Up to t 5 / 2 
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and to the third order of t: 

\^(x>L,t)\ 2 ~ -L x 

r i + 4 M +4 3( fco ,) 2 -5 f2+ A 8_ | 
[a; a; 4 ar x J 



we can see that the barrier's presence is felt only at 
the third order of t. 

Even when k ^ X and x — ► oo the barrier's presence 
has a significant influence when the measurement is 
taken in the range 4fcoa;/A 2 <iCi/A. 
On the other hand, for < x < L, owing to the 
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FIG. 2: A comparison between the solution with (lower panel) 
and without (upper panel) the barrier, which is represented by 
the horizontal line at x = 1. The initial state is represented 
by the dashed line. The parameters in this case are L = 1, 
k = 30, t = 0.04 and X = 3. 

reflection from the barrier, the probability density's 
dependence on the barrier appears even in the t 3 / 2 order. 



Moreover, the difference between the x > L and x < L 

regimes is clear from the figure, fn the latter regime the 

influence of the potential is felt for longer distances, but 
i i 2 

still when ^ > oo its influence decays to zero. 
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FIG. 3: : The difference between the solution when the barrier 
is present and when it is absent. Clearly, as \x\ 2 /t 3> 1 the 
two solutions are identical. 
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When the potential is absorptive the Schrodinger equa- 
tion should be rewritten 

d 2 dib 
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7T 
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That is, the dependence appears at the probability 
density at the coefficient of t 2 . 

Obviously, this approximation applies only when the 
argument L[L — x)/t is not too small. 

In Fig. 2 we plot a comparison between the propagation 
of the wavepacket in case the barrier is absent (upper 
panel) and when it is present (lower panel). In Fig. 3 the 
difference between the two (with and without the barrier, 
i.e., A|^| 2 = \ij}\l ith - \Mluhout) is Plotted. Despite 
the fact that the packet passes through the potential, 
its effect beyond the barrier x > L is miniscule and 
for \x\ 2 ^> 1 the two solutions are essentially identical. 
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and the probability density satisfies 



(12) 



\^(x>L,t)\ 2 = -t-\l + 2^^t + 0{t 2 )\. (13) 



In this case the potential presence appears in the proba- 
bility density at the second order of t (and not the third 
as in eqJSJ. 

In the case where the wavefunction initially vanishes at 

z = o (HQ3G3), e.g., when 



ip(x, t = 0) = 0(—x) sin(fcox) 
The general solution is 



(14) 



4 



a ix 2 /4t 



4? 



-W 



- ko) 



Ji r/it 



A/2 



Jx 2 /4t 



4l 



-w 



2t 



k - iX/2 

Jy 2 /4t X j 2 

~ k + i\/2 
I 



~ ko) 



2t 



'it[ i—) 

y 2t 2' 



2t 



(15) 



In short times one can expand this expressions in powers 
of t. 

Up to t 5 / 2 
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the conclusions. However, the main conclusions of the 
semi-infinite plane wave can be deduced even when the 
initial wavepacket is a continuous wavefunction, such 
as a Gaussian. However, owing to the finite spectral 
width of the Gaussian, high energy particles are very 
rare in the packet, and therefore the reasoning that was 
used in the semi-infinite plane wave can be used only 
in the intermediate temporal region and fails at t — > 0. 
Therefore, we should expect to find the same conclusions 
in a certain intermediate period (as in ref. ||J). 
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If the initial wavepacket is a Gaussian: 



and 



7T2; 4 



(18) 



we can see that the barrier is not felt even at the third 
order of t. 

Although the phase of ip , " feels" the barrier even in the 
first order of t, the leading term is proportional to 



then 



ifj{x, t = 0) = J-- exp |-(-) 2 + ik x) (20) 
V 7r a v a ) 
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III. THE CONTINUOUS CASE: GAUSSIAN 
DYNAMICS 



^(x>L,t>0)= (21) 

-i p exp i — (k — k ) 2 ct 2 /4 + ikx — ik 2 t\ 
dk ^ J 
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The fact that we used singular initial wavefunction 
may raise skepticism about the physical validity of 



the general solution is 
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where s = y cr 2 /4 + it. For short times 



tp(x > L,t > 0) 



> : aJ 5 exp 1 4 
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(23) 
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where for large distances x + L o 2 ko can be approxi- 



mated in the two extreme cases: t <§C cr 2 /4 and t ^> cr 2 /4. 
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In the former case 



tjj(x > L,t > 0) = 




the barrier influence is independent of time, and in fact, 
very far from the barrier x 3> Act 2 the barrier's influ- 
ence is negligible. However, we see here that in the limit 
t —>■ 0, due to the finite spectral width ~ a , the pres- 
ence of the barrier is always felt. (Note that the singular 
case 4x = Act 2 is not consistent with the above approx- 
imation). This is to be expected, since in a Gaussian 
distribution the number of particles with extremely large 
enegies is exponentially small. 
When t > cr 2 /4 

ip{x >L,t>0) = —== rrexp^ -2- \ 

we recognize a penetration velocity. When i/f< A the 
barrier has a large impact, however, if the particles' ve- 
locity is very large x/t ^> A the barrier's influence is 
negligible. 

And similarly, in the temporal period c 2 /4 -C t <C 
x/X the difference between real and imaginary barrier is 
apparent. For a real barrier 

,, 9 1- (Xt/x) 2 f er 2 fc 2 1 , s 
\^{x>L,t>Q)\ 2 = L^expj (25) 



while for an imaginary one 



.„ 1 - 2Xtlx { a 2 kl \ , . 
\4,(x>L,t>0)\ 2 = 2 J exp| (26) 



Again, the the influence of the absorptive potential 
appears in a lower order term. 



IV. SCHEMATIC EXPERIMENTAL 
REALIZATION 

One of the methods of emphasizing the impact of the 
potential is by placing the potential barrier (delta func- 
tion in our case) in one arm of a Mach-Zehnder interfer- 
ometer (see fig. 4). 

Let us denote by c\ and ic-i the transmission and reflec- 
tion coefficients of each of the two interferometers' beam 
splitters (BS1 and BS2 in figure. 4). With this notation 
we assume (without loss of generality) that both c\ and 
C2 are real and conservation of energy implies c 2 + c 2 , = 1 . 
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Let us further assume that both BS's are almost 50:50, 
i.e., c\ = 0.5(1 - e), cf, = 0.5(1 + e) and e < 1. In this 
case, for short time 

ip(x,t) ■ 



it e" 2 / 4i 
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and the probability density can be approximated 
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Thus, when the potential is real the potential-dependent 
term has a cubic dependence on time 



t 



\^{x,t)\ 2 ^ — 



X 
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(29) 



while if the potential is imaginary the temporal depen- 
dence of the potential-dependent term is parabolic 

|V>(x,t)| 2 ^^[e-— t] (30) 
irx z x 



To emphasize the difference we define A|?/>| 2 = l^j 2 — 
\ip\ 2 ree as the difference between the probability density 
at the interferometer exit when the barrier is present 
(\ip\ 2 ) and when it is absent (iV'l/ree)- 

In fig. 5 we plot the temporal evolution of A|^| 2 , which 
is measured at the exit of the interferometer (at x — 10 
from the barrier) for the two cases (real and imaginary 
potentials). The only difference between the two plots is 
the potential (iX instead of A). While the two plots are 
similar after long times, their temporal differences are 
considerably different for short times as eqs. 1291 and I3HI 
imply (like t 3 and t 2 respectively). 



V. SUMMARY 

The short-time influence of a delta-function potential 
barrier on an initially confined wave-packet was investi- 
gated. It was shown that at short times the barrier pres- 
ence has a negligible influence, if any, on the wavepacket 
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dynamics. This result applies also for the probability 
density of the particles that passed through the barrier. 
It was also demonstrated that at short times an absorp- 
tive barrier (i.e. imaginary) has a different impact on the 
dynamics than a non-absorptive (i.e., real) one. Namely, 
at short times an absorptive barrier appears at the co- 
efficient of the t 2 term, while a non absorptive barrier 
appears only at the coefficient of the t 3 term. Therefore, 
it is possible to distinguish between an imaginary and 
a real potential barrier by measuring its effect at short 
times only on the transmitting wavefunction. There is 
no need to measure the transmission and reflection coef- 
ficient simultaneously. 

FIG. 4: schematic illustration of a Mach-Zehnder interferom- 
eter for barrier classification. 
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FIG. 5: The temporal evolution of A\ip\ 2 outside the interfer- 
ometer (a distance x — 10 from the barrier) for a real potential 
(solid line) and an imaginary one (dashed line), cl = V0.49 
and the other system's parameters are the same as in fig .2 
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